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Abstract 

The role of Bose-Einstein correlations in a widely discussed inter- 
mittency phenomenon is reviewed. In particular, it is shown that par- 
ticle correlations of different origin are better displayed when analysed 
as functions of appropriately chosen variables. Correspondingly, if the 
shape of Bose-Einstein contribution is chosen to be Gaussian in 3- 
momentum transferred, it provides the power-like law in 4-momentum 
squared and is smeared out in (pseudo) rapidity. 

1 Introductory survey and definitions. 

The study of fluctuations and correlations in hadron production at high ener- 
gies has found considerable interest in recent years. The g-particle inclusive 
densities Pq{pi, ■ ■ ■ ,Pq) or rather the factorial moments (n(n— 1) . . . {n—q+1)) 
estimated in different phase space regions 6 were studied in a variety of re- 
actions ranging from e~^e~ to nucleus-nucleus scattering. 

The concept of intermittency has been introduced in order to describe 
enhanced fluctuations observed for individual events in the density distribu- 
tions of hadrons (for a recent review see [1]). Originally the definition of 
intermittency was the strict power laws of the normalized factorial moments 
taken as a function of bin size [2] at small bin sizes (self-similarity of the 
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moments). Later on it has become customary by calling intermittency any 
increase of the factorial moments with decreasing phase space intervals with- 
out regard to scaling behaviour. The extended version indicates just positive 
correlations in momentum space which increase at smaller bins. Such an in- 
crease of factorial moments demonstrates that the multiplicity distributions 
become wider in small phase bins i.e. fluctuations arc stronger. 

Being important by itself, it shadows a crucial feature of the original 
proposal, namely, the search for self-similarity in the processes which could 
be related [3] to the fractal structure of the pattern of particle locations in the 
three-dimensional phase space. Therefore using the term "intermittency" in 
a wide sense here, we still keep in mind that the scaling regime is of primary 
interest. It gave name to the whole effect when was first proposed in study 
of turbulence [4] . 

We shall discuss some possible sources of correlations. Our main concern 
is to show that they contribute differently when the various projections of 
multiparticle phase space arc considered. The proper choice of the variable 
can emphasize the particular correlation and, vice versa, conceal the contri- 
bution of others. Therefore, by choosing different variables we study the 
correlations of different origin weighted with different weights. 

Various explanations with different physical origin of correlations have 
been proposed. We shall discuss them in more detail in Section 3, but here 
let us just mention some of them. The analogy to turbulence has led in the 
original proposal to the term "intermittency" and to the phenomenological 
multiplicative cascade model of the phenomenon. On more strict grounds, it 
could be related to the parton shower in quantum chromodynamics (QCD) 
as described in Section 3.5. The observed phenomenon of hadron jets by 
itself indicates strong positive correlations. Thus, the scaling regime could 
appear also at the stage of transition from partons (quarks, gluons) to the 
observable particles (mostly hadrons) . Therefore the ideas of phase transition 
have been elaborated too. Even a more trivial dynamical reason could be 
connected with abundant production of resonances which, surely, imply the 
correlation of relative energies of final particles. Some unknown sources of 
completely new dynamics have been looked for (e.g. stochastic dynamics, 
instabilities etc.). 

Beside those "dynamical" effects there exists the well-known symmetry 
property of the interaction which necessarily contributes to the enhancement 
of correlations in small phase-space volumes. We mean Bose-Einstein corre- 
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lations due to the symmetrization of the wave functions of identical particles. 

In this review, we shall deal mostly with the problem of what and how 
the Bose-Einstein effect contributes to the observed increase of moments of 
multiplicity distributions in ever smaller phase space regions when various 
variables are chosen to tag the phase space bin. Other topics, mentioned 
above, will be considered to the extent they are necessary for clarification of 
some issues related to our main purpose. 

Bose-Einstein (BE) statistics leads to specific positive correlations of iden- 
tical bosons (in multiple production processes those are mainly like-charged 
pions). BE correlations are stronger for particles having smaller difference of 
three-dimensional momenta. So BE correlations lead to intermittency effect 
(in the broad sense of the word). Indeed it was first considered in [5], and 
was asserted a few years ago [6, 7, 8] that BE correlations may be responsible 
for the intermittency effect. However ineffective methods of the data analysis 
did not give a possibility to reach convincing conclusions that time. 

Early investigation of the intermittency was concentrated on the study of 
one-dimensional rapidity dependence of the traditional factorial moments in 
decreasing rapidity intervals. Then it was realized that intermittency effect 
in multiparticle production, if any, would naturally occur also in two and 
three dimensions when two or all three components of momenta of outgoing 
particles are registered. 

Furthermore high-order inclusive density Pq{pi, ■■■,Pq) contains contribu- 
tions from lower order correlations. Being interested in genuine multipar- 
ticle correlations one is led to investigate (connected) correlation functions 
Cq{pi, ...,Pq) with lower order correlations subtracted (the correlation func- 
tion Cq vanishes when a subgroup oim < q particles is statistically indepen- 
dent from the other q — m particles) . 

The normalized factorial moment estimated in the phase space region 5 
is defined as 

Fq{S) = = ^J/Pl...l dpqpq{p,^ . . . , Pq) ^ 

_ (n(n-l)...(n-g + l)) 

where F^'^^ are called unnormalized moments. The normalization is chosen 
so that at integer ranks q one gets Fq = 1 for Poisson distribution. In Eq. (1) 
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dpj is the differential momentum space volume 

dpj = d^pj/2Ej{2nf (2) 

or any other differential interval of interest. Here 5 will be used as a general 
notation for a selected phase space vohime which differs for different selec- 
tion procedures. In one-dimensional analysis this is usually the interval of 
rapidity 5y or pseudorapidity 5r] where rapidity is defined as y = | In 
and pseudorapidity as 77 = |ln|^ = In tan | with E,p,pi,9 being the en- 
ergy, the momentum, the longitucfinal momentum and the angle of particle 
emission. 

In practice, an averaging over cells of the original phase space is per- 
formed, so that traditional normalized factorial moments in d dimensions are 
determined as 

Pg^^) = ^dT.Km{S). rf= 1,2,3, (3) 

m=l 

where the sum runs over M"^ dg-dimensional boxes having the same size. 
Self-similarity of the moments Fq{5) taken as a function of momentum cell 
size suggests a power law behaviour 

^.(^)' (4) 

where (pg > are known as intermittency indices. 

General interrelation of the inclusive densities pq and the correlation func- 
tions Cq is provided through the inclusive generating functional 

00 1 <? „ 

G[Z] = E n / dPjPM, ■ ■ ■ ,Pg)Z{pi) . . . Z{pq), (5) 

q=0 «■ j=l •' 

where Z{pj) are auxiliary functions. The inclusive densities are then given 
by differentiation 

'^^^'-■'P^^ = 6z{p,y^.szipq)^'-'- 

The correlation functions Cq are defined through the generating functional 
G in the following way (cluster expansion): 

G[Z] = exp( / dp,p^{p)Z{p) + E -T n / dpjCqip,, . . .,Pq)Z{p,) . . . Z{pq)) 

(7) 
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being "the exponents" of the inclusive densities, 

5" In G , 

A comparison of (5) and (7) leads to relationships 

P2{PUP2) = Pl{Pl)pl{P2) + C2{pi,P2), (9) 

P3{Pl,P2,P3) = Pl(Pl)Pl(P2)pi(P3) + 

3 

J2 Pl{Pi)C2{Pj,Pk) +C3{pi,p2,P3) (10) 

etc, or the inverse ones: 

C2{pi,P2) = p2{pi,P2) - Pi{pi)pi{p2), (11) 

3 

C3{PUP2,P3) = P3{PI,P2,P3) - Pl{Pt)p2{Pj,Pk) 

+2pi(pi)pi(p2)pi(p3) (12) 

etc. Let us note that the correlation functions have the above mentioned 
physical interpretation if the number of produced particles exceeds (not clear 
a priori to what extent) the order of the correlation function. 

Integrated (and normalized) correlation functions are known as cumulant 
moments 

ii:(")(5) 1 /• r 

Cumulants of Poisson distribution are identically equal to zero. 

The generating function of the moments arises if one puts the functions 
Z{pj) in generating functional (5) to be a constant z, 

oo -| oo 1 oo 

9{^) = E = exp((n)z + Y: -K^^^z'^) = ^n(l + ^)" (14) 

g=0 5=2 9- n=o 

The moments and particle number distribution P„ can be found by differen- 
tiation of the generating function: 

^^|.=o = i^i"^ = {nYF, = {n{n - 1) . . . (n - g + 1)), (15) 



*'"»<'-\.„ = if<«) = („>.if„ (16) 
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,.-. = P.nl (17) 

Relationships between factorial moments and cumulant moments follow from 
(14) or from (10) and (12): 

^2 = 1+^2, (18) 

Fs = 1 + 37^2 + Ks, (19) 
F4 = 1 + 6i^2 + ^Kl + 4Ks + i^4, ■ ■ ■ (20) 

for a fixed cell in momentum space. An averaging over cells, such as that in 
(3), 

K<l = W^i: ^^^rniS) (21) 
m=l 

requires the corresponding averaging of Eqs. (20). 

Let us return to the intermittency study. The cumulant moments (21), 
representing genuine multiparticle correlations, suffer from arbitrary binning 
and low statistics. The accuracy in the measurement of the factorial moments 
is also unsatisfactory for small bin sizes. That has led to investigation of more 
general density and correlation integrals which give the possibility to use the 
available statistics in an optimal way and to perform the correlation analysis 
in any convenient variable. 

The general density integral is defined as 

F^i^) = F n / dPjPM, ■ ■ ■ ,pMS;pi, . . .,P,), (22) 

Q j=l 

and the correlation integral as 

^? = 1^ n / dPjCgiPi, . . .,Pg)n,{5;p,, . . . ,pg), (23) 

^^<7 j = l ■> 

where the "window function" Vtq is determined to be nonzero in some pre- 
scribed interval of the g-particle phase space. The normalization to uncorre- 
lated background is suggested in (22), (23): 

^1=11 I dpjpj{pj)ng{5;pi, ...,Pj). (24) 
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Table 1: Intermittency indices for all charged particles and negatives only 
(for small bins ) 





02 


03 


04 


05 


all charged 


0.008 ± 0.002 


0.043 ± 0.006 


0.16 ±0.02 


0.39 ±0.06 


negatives 
only 


0.007 ±0.003 


0.06 ±0.02 


0.29 ±0.06 





The general definitions (22), (23) give a possibility not to split the g-particle 
phase space in an artificial way (as it was the case in Eqs. (3), (21)) ensuring 
better statistics. This in turn gives a possibility to measure and for 
like-charged and unlike-charged particles separately providing direct evidence 
for the BE correlations. 

2 Experimental survey. 

The extensive review of experimental data is given in [1]. Here, we mention 
some typical findings which help us exemplify our treatment, especially, those 
which reveal BE-contributions. 

2.1 One-dimensional moments (rapidity variable). 

Consider first the one-dimensional analysis in rapidity variable. Factorial 
moments (3) for like-charged particles were measured in 7r~^p and K^p inter- 
actions at 250 GeV/c by NA22 Collaboration [9]. They are characterized by 
intermittency indices 0q according to the fit 

F,{6y) = a.iSyy'^^. (25) 

Though the bin size dependence turns out to be far from power-like in the 
whole interval measured, but it could be fitted by power laws separately at 
large {S > 1) and small {S < 1) bins. We remind that the latter region was 
suspected for such a law in original proposal (see Table 1 and Fig. 1). 

The values of the moments are smaller but the increase with decreasing 
bin sizes is stronger when a like-charge sample is used instead of all charged 
sample. To what extent the BE correlations (presumably dominating the 
correlation of like-sign particles) are responsible for the full intermittency 
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is not clear from the data in the rapidity variable. Anyhow they are not 
dominating. 

Analogous analysis of the factorial moments was performed in pp colli- 
sions at = 630 GeV by UAl Collaboration [10], sec Fig. 2. The authors 
conclude that the BE effect is weak in this variable. (Their early data of 
intcrmittcncy were analysed by moans of the negative binomial distibution 
and the pure birth stochastic theory in Rcf.[ll].) 



2.2 Higher-dimensional analysis of the factorial mo- 
ments. 

The intermittency effect is more pronounced in two and especially in three 
dimensions. The impressive results on the second factorial moment were re- 
ported in i^N interactions at 280 GeV/c by NA9 Collaboration [12]. The 
moments F2 for unlike- charged and for negative particles were measured in 
one and three dimensions. The number of boxes is M and correspond- 
ingly. The results are given in Fig. 3. The authors claim that the strong 
intermittency signal observed in F2 in three dimensions has to be attributed 
exclusively to BE correlations since such a signal is not present in . This 
conclusion is supported by the fact that the Lund model (not containing BE 
correlations) is in rough agreement with the data for but in complete dis- 
agreement for F2 . In one (rapidity) dimension the effects are much weaker. 
Nevertheless a considerable difference between data and the Lund prediction 
is seen for F2 in one dimension too whereas for F^~ the difference is small. 
This also supports the noticeable role of BE correlations. 

Similar results for F2{d) in three dimensions were also presented in 7r~^p 
and K^p interactions by NA22 Collaboration [9]. The box volume depen- 
dence was fitted with the form 

K.^{d) =F2{6)-l = c + a5-\ 5 = (26) 

A striking difference for unlike- and like-charged pairs was found. While 
the ( H — ) pairs are dominated by long-range correlations (large c), these 

are smaller or absent in the case ( ) and (-1- -|-). Correspondingly, the 

parameter a is compatible with zero for (-1 — ), but relatively large for ( — ) 
and (H — h). This again supports BE interpretation of the intermittency in 
three dimensions. 
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2.3 Density integrals. Three dimensional analysis. 

For higher orders (g > 3) the traditional normahzed factorial moments have 
large statistical errors when evaluated for small bin sizes. So nowadays the 
general correlation analysis is usually performed (see Eqs. (22)-(24)) per- 
mitting not to split the phase space interval under consideration. Using 
the density and correlation integral method, one must introduce the dis- 
tance 5ij for each pair of particles and confine the permissible g-particle 
phase space. Usually the "distance" is defined as four-momentum difference 
squared, 5ij = Q^j = {pi —pjY- The permissible phase space is most com- 
monly confined either by GHP condition [13], Sij < S for all i,j<q or by so 
called star topology, 5ij < S ior j — 2, . . . ,q, or by snake topology, Sjj+i < S 
for j = 1, . . . , q — 1. 

A comparison of GHP-integral with the conventional normalized factorial 
moments was performed in vr+p and K~^p interactions at 250 GeV/c by NA22 
Collaboration. The three-dimensional "distance" between two particles i and 
j was defined using y, (p and pt variables in a rather specific way as 

dij = max{\yi - yj\, \ipi - ipj\, \pt,i - Pt,j\)^. (27) 

The size S of g-tuple was defined by the smallest box volume that encloses 
the g-tuple. The determination of the density integrals can now be compared 
with moving around a box in the full phase space under consideration and 
counting the number of the g-tuples fitted into the box. 

In Fig. 4 conventional moments F2 and F3 are compared to the density 
integral version F^^^^ and F^^^ . As anticipated, in Fig. 4 one indeed 
observes that statistical errors in the F^^^ are strongly reduced. This, in 
principle, allows the analysis to be carried down in much smaller box volumes. 
It, furthermore, allows a comparison of different charge combinations. 

For the second-order integral F^^^{5) the fit (26) was used as it was the 
case for conventional factorial moments considered in the previous section. 
The results are shown in Table 2 (we omitted the constant b when the con- 
stant a was found to be compatible with zero) being qualitatively the same 
as for ordinary three-dimensional factorial moments. 

The above results again support the conclusion that the intermittency in 
three dimensions is strongly enhanced due to BE correlations. 

For the higher order GHP density integrals the modified power law as- 
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Table 2: Results of fits (26) to the data on F^"^ 





unlike charged 


negatives only 


positives only 


a 


0.0006 ± 0.0009 


0.0115 ±0.0003 


0.04 ±0.02 


b 




0.469 ± 0.004 


0.37 ±0.06 


c 


0.380 ± 0.006 





0.08 ±0.03 



Table 3: The slopes aq /a2 obtained by fitting (28) to the data 





all charged 


negatives only 


positives only 


a3/tt2 


3.81 ±0.09 


4.3 ±0.2 


5.3 ±0.2 




8.0 ± 0.03 


9.4 ± 0.5 


15.0 ±0.5 


a^loL2 


11.8 ±0.3 


13 ±0.1 


24 ± 1 



sumption 

lnFf^^((5)=a,±^F2((5) (28) 

OL2 

can be fitted to the data though (28) only holds approximately. The results 
are presented in Table 3. 

As can be seen from the Table, the growth in decreasing three-dimensional 
phase-space volumes is faster for higher order density integrals and for like- 
charged particles. 

Let us mention here the analysis [14] done by the same NA22 Collabora- 
tion using the opening angle Qij between two particles 

% = arccos(piPj/|pi||pj|) (29) 

with pj and Pj being the three- momenta of particles i and j. An angular dis- 
tance measure for more than two particles is defined as the maximal relative 
angle between all the pairs chosen. Therefore, the numerator of the factorial 
moment of rank q is determined by counting, for each event, the number of 
g-tuples that have a pairwise angular opening smaller than a given value and 
then averaging over all events. 
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The fitted values of intermittency indices appeared to be very low (some- 
times negative) and strongly dependent on the production angle of the parti- 
cles. We argue in Section 3.5 that such an analysis reminds of two-dimensionall 
analysis in relative pseudorapidity and azimuthal angle but with the bin size 
which depends on the production angle. Therefore it is not very useful for 
comparison with analytical calculations for parton jets. Nevertheless, it is 
possible to compare the experimental findings with the FRITIOF Monte- 
Carlo model. It shows that BE correlations should be incorporated into the 
model to get the agreement. 

2.4 (5^-analysis of the density integrals. 

NA22 Collaboration has also presented the data on Q^-dependent GHP den- 
sity integrals, that is the integration over particle densities pq was confined 
by the factor 

^^."""m-fim'-Qi)- (30) 

The data are shown in Fig. 5 and Table 4 where the parameters of the power 
law fit 

F^'^'^iQ') = a,{Q')-'^^ (31) 

are given. One can see that evaluation of the GHP integral for different 
charge combinations yields effective intermittency indices a factor 1.2 (^5) 
to 1.6 (^2) larger for the negatives-only sample than for all-charged sample. 
This indicates the important role of BE correlations as displayed on Q^-scale 
though more definitive conclusions can be hardly reached from the data. 

Let us note in this connection that pj-dependence of the second and third 
order GHP density integrals taken as a function of Q"^ appears to be opposite 
for all-charged and negative-only particles in this experiment, see Table 5. 

According to the Table, the intermittency effect is weaker at low pt than at 
high Pt for negatives. On the contrary, intermittency is stronger when small 
Pt particles are selected from all charged. This hampers easy interpretation 
of the all-charged data only, without proper analysis of the different charge 
combinations. 

The related results on density integrals F2{Q'^) and F^,{Q'^) were presented 
in fiN interaction at 490 GeV/c, using data from the E665 experiment at the 
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Table 4: Results of fits to the data presented in Fig.5 according to Eq.(31). 





all charged 


negatives only 


positives only 


02 


1.219 ±0.003 
0.051 ±0.001 


1.131 ±0.002 
0.081 ± 0.001 


1.026 ±0.002 
0.067 ±0.001 


as 

03 


1.751 ±0.007 
0.177 ±0.002 


1.38 ±0.01 
0.253 ±0.004 


1.15 ±0.05 
0.227 ±0.003 


04 


2.90 ±0.02 
0.358 ±0.006 


1.88 ±0.02 
0.47 ±0.01 


1.41 ±0.01 
0.45 ±0.01 


05 
05 


5.45 ± 0.08 
0.56 ±0.01 


2.78 ±0.07 
0.66 ±0.02 


1.89 ±0.04 
0.66 ±0.02 



Table 5: Intermittency indices for different pt regions 







Pt < 0.15 GeV/c 


Pt > 0.15 GeV/c 


all-charged 


02 


0.046 ± 0.002 


0.032 ± 0.001 




03 


0.136 ±0.004 


0.107 ±0.003 


negatives-only 


02 


0.053 ±0.002 


0.081 ±0.001 




03 


0.17 ±0.01 


0.269 ±0.006 
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Tevatron of Fermilab [15]. This time the star topology was used to confine 
the g-particle phase space, that is the confining factor was taken in the form 

^g(Q')-fim'-Qij)- (32) 

The Fig. 6 shows a log-log plot of the density integrals for different charge 
combinations. F2{Q'^) rises more steeply with for ( — ) than for ( H — ) 

pairs. The same is true for F^^Q"^) for ( ) triplets as compared to (ccc) 

triplets. This different behaviour indicates important role of BE correlations 
between hke-sign particles. 

No significant energy ll^-dcpendence of the F2 slope was observed for 

> 0.01 GeV^ ; for smaller the F2 slope of ( ) pairs seems to be 

somewhat larger for the high-iy sample. In fact, the slopes (ilnF2/(iln(l/(5^) 
of the F2 integrals show close agreement for NA9 and E665 experiments 
[12, 15] in spite of the somewhat different (1^) values and the experimental 
differences of the two experiments. 

Q^-analysis of the density integrals was also performed in pp reaction at 
§1/2 = 530 GeV by UAl Collaboration [10]. The authors however used quite 
another definition to confine the permissible g-particle phase space, 

t < (33) 

ii<i2 

proposed in Ref. [16], which is much stronger than Eqs.(30),(32). The results 
are given in Fig. 7 for all-charged and like-sign particles up to order q=b. 
The Q^-dependence of the density integrals is close to linear one in a log-log 
plot. The corresponding intermittency indices are listed in the Table 6. 

The comparison in Fig. 7 and Table 6 shows once again that the slope 
parameters (pq are bigger for like-sign particles than for all charged particles 
and the condition 

Cpqilike — sign) ~ 2(pq{alV) (34) 

is fulfilled approximately in the (^^-representation with the confinement con- 
dition (33). Let us note, that the relationship (34) was suggested in Ref. 
[7] as an indication on BE origin of the intermittency effect (though it was 
advocated for rapidity variable in [7]). 

At the same time only small differences between different charge com- 
binations are visible in the pseudorapidity analysis of the same UAl data. 
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Table 6: The results of fitting of normalized density integrals as functions 
of to a power law in pp interactions at s^/^ = 630 GeV 



slope 

parameters 


02 


03 


04 


05 


all charged 
particles 


0.0348 ± 0.0006 


0.078 ±0.001 


0.213 ±0.004 


0.338 ±0.019 


like-sign 
particles 


0.0522 ± 0.0009 


0.147 ±0.001 


0.443 ± 0.01 


0.855 ± 0.051 



see Fig. 2. The Figures 2 and 7 demonstrate that the manifestation of BE 
correlations is strongly dependent on the variable used. 

2.5 Second order correlation function. 

A detailed study of the second order normalized correlation functions C2(<5^) 
(sometimes known as differential correlation integral) was undertaken in 
interaction by UAl Collaboration [10]. Fig. 8b shows a comparison of the 
samples of the like-sign pairs with unlike-sign pairs and all charged particles. 
In Q^-variable one observes a dominance of unlike-sign pair correlation for 
0.03 GeV^ < <1 GeV^ which is at least partly due to resonances and 
particle decay (i.e. there is a broad peak at ~0.5 GeV^ which is due to 
p- meson decays and a peak at Ri0.17 GeV^ which is due to K\ decays). 
However at small <0.03 GeV^ this function is nearly constant. Contrarily, 
the like-sign particle correlation function is rising very rapidly at very small 

up to = 0.001 GeV^ exhibiting approximately power-law behaviour. 

A comparison with the same analysis in pseudorapidity variable 6r] (Fig. 
8a) shows once more the significance of the choice of the proper variable in 
correlation analysis: the two body correlation function of all charged particles 
is dominated by unlike-sign particle correlations when analysed in 5r] but 
dominated by the like-sign correlation function when analysed at very small 
Q^. Similar results on C2{Q'^) were reported in vr+p and K^p interactions by 
NA22 Collaboration. 

The correlation function C2{Q'^) was also investigated in e'^e~-annihilation| 
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at 91 GeV (Z°-boson) by DELPHI Collaboration (see [17]). In the range 
g2 <0.03 GeV^ the function C2{Q^) of DELPHI (e+e") and that of UAl 
(pp) show quite a similar shape, see Fig. 9. Since there is a rise for smaller 
only for the same charged pairs, BE correlation is evidently responsible 
for this behaviour. 

However, the C2 of DELPHI is also rising for smaller in the interval 
> 0.03 GeV^ both for same and opposite pairs, where the UAl data show 
a comparatively small rise. Therefore in e'^e~ annihilation some other mech- 
anism must be responsible for this power law behaviour which is manifest 
even in oppositely charged pairs. Jet evolution or hadronization may play 
the role. 

2.6 Correlation integrals. 

A further insight into the problem of correlations would be provided by in- 
vestigation of cumulant moments Kg{6) describing true multiparticle corre- 
lations of the order q with lower order contributions subtracted (cumulants 
vanish whenever one of the particles involved is statistically independent of 
the others). The cumulants relevant to intcrmittcncy study are integrals over 
(connected) correlation functions taken in decreasing phase space volumes. 
An experimental investigation of the correlation integrals is difficult because 
they suffer from low statistics. It was performed recently using star topology 
confinement in correlation integrals and g?- as a "distance" between particles, 
that is inserting the factor (32) into the correlation integrals. 

Results on correlation integral K^{Q'^) for all charged particles and neg- 
atives only in interactions [15] are shown in Fig. 10a. The behaviour 
of K^iQ"^) is approximately power- like in variable. In order to find the 
origin of the three-particle correlations in Fig. 10a, jiN Monte Carlo (MC) 
events were generated according to the Lund model without BE correlations. 

The MC predictions for K-^^Q"^) are shown in Fig. 10b. For ( ) triplets, 

K^iQ'^) of the MC events is rather independent of l/Q"^ in contrast to the 
data. This shows that the rise of K^iQ"^) in the data is very likely due to 
three-particle BE correlations which were not incorporated into the Lund 
MC used. For (ccc)-triplets the situation is more complicated since in the 

data both BE correlations (weaker than in ( )) and resonance decays 

contribute. In the MC (without BE but with resonance decays) K^iQ'^) is 
smaller than in the data but rises due to resonance decays. 
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Table 7: Slopes (f)q fitted to the dependence Kq ~ (Q^) 



charge 
combination 


q = 2 


q = 3 


q = 4 


q = 5 


all 


0.205 ± 0.005 


0.72 ± 0.03 


1.2 ±0.2 


2.0 ± 1.0 


like-charged 

++ 


0.387 ±0.009 

0.438 ±0.010 


1.03 ±0.08 


1.8 ±0.3 




unlike- 
charged 


0.096 ±0.004 


0.61 ±0.03 


2.0 ±0.5 





Irreducible higher-order correlations were also established up to fifth or- 
der in multiparticle production in 7r+p and K^p collisions at 250 GeV/c by 
EHS/NA22 Collaboration [18]. The star integral method has provided a clear 
improvement over the earlier analysis based on the same data. The charge 
dependence of the correlation was studied in a comparison of like-charged 
(Fig. 11a) and unlike-charged (Fig. lib) particle combinations. Both charge 
combinations show non-zero genuine higher order correlations and an increase 
of the correlation functions with decreasing interval . 

The correlations among unlike-charged combinations (i.e. combinations 
to which resonances contribute) are relatively strong near ~1 GeV^ but 
the increase for smaller is relatively slow. Correlations among like-charged 
particles are small at ~1 GcV^ but increase rapidly to reach, or even cross, 
those of the unlike-charged combinations at lower . This difference dimin- 
ishes with increasing order q. The effective slopes of the power- like scaling law 
for the various charge combinations fitted in the range < ln(lGeV^/ Q'^) < 5 
are given in the Table 7. 

One can sec that like-charged particles show faster growth of the correla- 
tion integrals with decreasing Q^-intcrvals thus indicating the important role 
of irreducible higher order BE correlations. 

2.7 Summary of experiment. 

Correlations between different charge combinations in multiple production 
processes were recently measured in a variety of reactions {pp, nN,e^e^ , 
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n^p, K^p). The data on two- and many-particle correlations taken as a 
function of permissible phase space volume were presented (intermittency 
study). This new development became possible due to investigation of the 
density and correlation integrals which give a possibility to use the available 
statistics in an optimal way and introduce any convenient variable. 

From this analysis it becomes clear that the increase of correlations with 
decreasing 3-dimensional phase space volume is essentially due to correla- 
tions between like-charged particles. An evident candidate for these like-sign 
correlations is the BE statistics. 

The intermittency effect (in its wide sense) depends strongly on variable 
chosen. It is poorly seen in rapidity-analysis and much more pronounced 
in three-dimensional- and Q^-analyses. The last variable is one of the most 
popular nowadays but it brings its own problems as will be discussed below. 
As we show, the effect in higher dimensions becomes stronger just due to the 
trivial decrease of the available phase space (smaller denominator in normal- 
ized factorial moments) while the steep rise at smaller is determined by 
the peaked contribution of BE correlations as exposed in that variable. 

3 Correlation studies in different variables. 
3.1 Variables and windows. 

For a proper interpretation of the intensity and correlation integrals and 
one has to look more carefully on their structure. In particular it is 
important to realize what kind of windows they provide. 

In general, the density and correlation functions Pq{pi, ■ ■ ■ ,Pq) and 
Cq{pi, . . . ,Pq) depend on 3g independent variables, say, 3(g— 1) independent 
momentum differences q^j = Pi — Pj and three components of their average 
momentum p = ^X^^^^pj. The 3{q — 1) differences must be confined in 
some or other way to provide an elementary cell in the phase space (the win- 
dow) and the remaining 3 variables show the position of the window. The 
dependence on the cell size is the object of the intermittency study. The 
dependence on the cell position is asubject of an averaging. The last is nec- 
essary to get an acceptable statistics. This averaging is explicit in traditional 
(vertical) factorial moments (3) and implicit in the density and correlation 
integrals (22), (23). Let us note that the cell size must be independent of the 
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order of the moments; otherwise they will not be the (averaged) moments of 
any distribution. 

The choice of appropriate variables is of great importance in correlation 
study. As it could be seen from the experimental survey the intermittency 
effect is poorly seen in rapidity variable and it is clearly seen in 3 dimensions 
and in Q^-variable. In BE correlations one is led to consider the differences 
of three-momenta of the particles qjj = Pi — Pj as input variables because 
the strength of BE effect is determined just by these differences. Considering 
other variables one has to translate the BE correlation from to these 
variables to compare the effect with the data. 



3.2 Q^-confinement. 

Let us consider in more details the Q^-variable which is intensively used in 
the experimental study of the intermittency. The "distance" between two 
particles is defined now as 

4 = (P^ - P^)' - i^i - Eif = - (35) 

and the factor confining permissible phase-space volume contains the corre- 
sponding (step-) ©-function; ©(Q^ — q^j)- It is not difficult to see that the 
above ©-function confines an ellipsoid 

+ - Q"^ (36) 

where M is the invariant mass, 

= g2 + 4^2^ and E is the total energy 
of the pair. So the confined volume in momentum space is 

The longitudinal phase space (along the direction of the total momentum 
Pi + Pi) is rising here with E and Q, 

= (38) 

On the one hand the ii^-dependence of V means that the Q^-confined cells 
(bins) have different sizes in momentum space being much larger for fast 
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particles. So with fixed one already has an averaging over momentum 
cell sizes. If the correlation functions have the momentum difference as a 
characteristic scale (as it is the case for BE correlations) then the above 
feature is rather unpleasant because the intermittency study suggests an 
investigation of correlations taken as a function of the cell size. 

On the other hand the particles in increased momentum intervals may 
well become uncorrelated and the integral over their correlation function is 
already saturated at some fixed momentum scale whereas the cell size Li is 
still rising with Q due to large E/M factor. In this case the normalized cor- 
relation integral will decrease with increasing (increase with decreasing 
Q^) approximately according to the power law in some interval. Just this 
mechanism was found responsible [19] for the steep rise of the two-particle 
correlation function of the like-charged particles with decreasing Q^. As a 
conclusion, variable is not an appropriate variable for the physical inter- 
pretation of the intermittency effect if the correlation between particles has 
a characteristic scale in momentum difference of- as it is the case for BE 
correlations (at the same time being related to invariant mass squared, 
Q2 _ ]\/[2 _ ^uil., is an appropriate variable for correlations arising due to 
resonance decay). The steep and approximately power law rise of the nor- 
malized correlation integrals with decreasing does not necessarily reflect 
the corresponding behaviour of the original correlation function Cq{pi, ...,Pq), 
being a kinematical effect inherent in variable. 

3.3 Three-momentum difference confinement. 

In general, different mechanisms responsible for particle correlation bring 
their own natural variables (the variable is "natural'' if the correlation func- 
tion Cq has a scale parameter related to this variable). If one considers BE 
correlations then the natural variable is three-momentum difference q^j = 
Pi "~ Pi (the energy difference is also involved for nonstationary particle 
sources and a modification of the variables is necessary for expanding sources).] 
The correlation function Cq{pi, . . . ,Pq) depends also on total momentum 
P = Z]j=i Pj (or on average momentum p — P/q). 

The essential point is that these two dependences have different momen- 
tum scales: the characteristic momentum difference scale p^ ^ 1/R {R is 
an effective source size for BE correlations) is noticeably smaller than the 
characteristic transverse and longitudinal momentum scales pt and pi. This 
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means that the correlation function is nonzero in a "tube" in 3g-dimensional 
space with its axis directed along the line pi = . . . = pg and with cross- 
section of the tube of the order The length of the tube is of the order 
Pi in longitudinal (along the beam direction) and pt in transverse direction 
{pi ^ Pt in experiment). 

Let us now confine the 3-momentum differences imposing the condition 
\^.^j\ ^ ^ (the window). The result for the normalized correlation integral Kg 
(see (23)) depends strongly on the relationship of S and the scale parameters 
Pd<.Pt <.Pi- 

a) The region of very small 5 < pd, where the correlation function is 
maximal, is not still accessible in experiments. 

b) In the region Pa < 5 < pt the integration over momentum differences in 
Kq is saturated with the scale Pd- All corresponding integrations are cut by 
5 in the normalization factor Nq. As a result, we are left with a dependence 
of the form 

c) In the region pt < d < pi only integrations over longitudinal momenta 
are cut off and we have a dependence 

d) In the full phase-space {6 ~ pi) the normahzed correlation integral 
takes its minimal value ^ 

^.~f#y • (41) 

The above rough estimation shows that the normalized correlation integral 
obeys approximately power-law 5-dcpendence in some subintervals of the 
window width 6. The dependence is steeper for small windows (the effective 
intermittency index in the region b) is three times larger than in the region 
c)) and for higher order of the moments. 

We conclude that the intermittency effect (in the broad sense of the 
word) is strong and ensures approximately power-like behaviour in three- 
dimensional analysis in momentum difference variables. No specific physical 
phenomena were involved to get the above power-law behaviour. It is just 
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connected with the trivial dependence of moments on the available phase 
space and does not involve any dynamical "anomalous'' dimension as it is 
the case for QCD jets discussed in the next section. The transition from 
"g^-regime" in (39) to "g- regime" in (40) is also due to the transverse mo- 
mentum limitation which is an inherent property of multiparticle production. 
Actually, the dependence on S is mostly determined by the "kinematical" 
normalization factor in the denominator which shows how large is the aver- 
age multiplicity within the phase space window provided by the confinement 
condition. 

3.4 An illustrative example (BE correlations). 

To get somewhat more detailed picture of the intermittency effect induced 
by BE correlations let us consider the cumulant moments for BE correlations 
in a simple environment. We suggest that the particles are created by some 
random gaussian currents [20]. ( A generating functional for BE correlations 
with chaotic and coherent components is formulated in Ref.[21].) Then the 
particle densities pq{pi, . . . ,pq) and the correlation functions Cq(pi, . . . ,pq) 
can be expressed through a single quantity F{pi,pj) which is an averaged 
current correlator. In particular 

Pi{p)^F{p,p), (42) 

and the (unnormalized) correlation integral (23) is 

^J"n<^) = (?-l)!n / dpjF{p,,p2)F{p,,p,)...F{pq,p^)Qq{S;p,,...,Pq). 

(43) 

We represent F-function in the form 

F{pi,P2) = [Pi{pi)pi{p2)]'^'du{pi-P2), d{0) = 1 (44) 

suggesting that the function dij depends only on the momentum difference. 
Remembering that the correlation scale pd inherent in dij is much smaller 
than characteristic particle momenta pt,Pi, we take the densities pi{pj) in 
Eq. (43) in coinciding points to get 

K^''\S)^iq-l)\ / dpp'i{p)lldipj-pj+i)dr2...dq,nq{5;p^,...,pq) (45) 
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with q — 1-fold integration over differences pj — Pj+i- 

Being interested in qualitative results we may cancel one of d-functions 
in Eq. (45) taking dqi — dqi(O) = 1. This leads to inessential numerical 
misrepresentation of the integral retaining its qualitative behaviour because 
one of q rf-functions in (45) does not serve as a direct cut factor. On an equal 
footing wc may neglect a variation of one-particle density pi {p) in the region 
where the density is substantial, this region being confined by pt and pi as 
above. The momentum difference window may be defined by a confinement 
of successive momentum differences, 

0,(5) = Q{5 - qi2)e{S - 523) ...e{S- g,-i,,) (46) 

(the snake topology). As a result, the correlation integral (45) takes a simple 
form 

/{pt iPi ) r 
dppl{J^dqd{q)y-\ (47) 

The normalization factor Ng in Eq. (23) taken in the same approximation is 
equivalent to (n)*, 

ript,pi) ^ fipt,Pi) 
NqPi j dPPiiJg dqy-^p^{n)\ (48) 

and the correlation integral (47) can be written in the form 

K'^^\5)^{q-l)\{nYk'-'' (49) 

with 

hdpd{p) jmin(prf,5) - ■ y ) 

The cumulant moments (49) correspond to generating function (14) of the 
form 

ln^(^) = -Hn(l-^). (51) 
This is the generating function of the negative binomial distribution (NBD) 

r(n + i) + + ■ ^^^^ 
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(let us note that the hnked pair approximation also leads to NBD [22]). The 
parameter "/c" of NBD is given by Eq. (50). It decreases with decreasing of 
permissible interval 5. In a rough step function approximation it reduces to 



in an accordance with estimations (39)- (41) of the normalized correlation 
integral. 

From the above example one can clearly see the role of BE-correlations in 
the intermittency effect. For very narrow windows (5 — > 0, A; — > 1) we have a 
geometrical distribution (GD) characteristic for thermal ("totally chaotic") 
excitation of the particle source. This distribution, having normalized cu- 
mulant moments Kg = (g — 1)! and normalized factorial moments Fg = q\, 
is the widest (having maximal moments) distribution possible in the present 
scheme, where BE correlations were estimated for a single fixed scale particle 
source. To get larger values of moments (for experimental evidence see Fig. 
3, Ref. [12]) one has to introduce additional fluctuations or to revise the 
particle source form. 

For wide enough windows (5 — > pi) one has A; ^ 1 leading to the relatively 
narrow Poisson distribution {PD — liuik^^NBD). This means that BE 
correlations which we consider here are not effective in large phase space 
volume. The change of the window width from zero to infinity interpolates 
between GD and PD and this is a general feature of the present scheme. NBD 
is one, especially simple, kind of the interpolation and it is not surprising that 
NBD appears as a possible approximation. More accurate estimations of the 
correlation integral (43) lead to similar conclusions. Additional coefficients 
7g appear in the cumulant moments (49) and a single parameter k varies 
slightly for different orders q, but the qualitative behaviour (49), (50) remains 
unchanged. 

3.5 The (pseudo) rapidity confinement. 

Let us turn now to the rapidity analysis which initiated the whole story of 
intermittency in particle physics. We show that pseudorapidity (which coin- 



k = l for 6 < pd, 
k S^/pI for pd<S < pt, 
k pa PtS/pl for Pt < 6 < pi, 
k ^ PtPi/Pd for 6 r^pi 



(53) 
(54) 
(55) 
(56) 
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cides with rapidity for rclativistic particles) is the most suitable variable to 
study correlations providing jet-like structure of high energy events. More- 
over, such events give rise to (quasi)intermittent (in strict sense!) behaviour 
of correlations in this variable. 

Really, we will speak about the one-dimensional angular confinement 
when the bin of a definite size 59 in a polar angle is chosen. The pseu- 
dorapidity size 5ri is proportional to 56 

5 = 5r]!^ 56/6 (57) 

with a weight given by a location 6 of the center of the bin that is trivially 
accounted when averaging over all locations is done. The dimension of the 
analyzed bin enters the results in a very simple (even trivial) manner as it is 
shown below. 

Let us consider in QCD the correlation within the gluon jet emitted by 
a quark produced in e"'"e~-annihilation [23, 24, 25, 26]. We are interested in 
correlations among the partons (mostly gluons) created during the evolution 
of the initial gluon and fitting the pseudorapidity bin 5r]. They belong to 
some subjet inside the primary jet which is separated from others due to 
the angular ordering in QCD. Surely, the prehistory of a jet as a whole is 
important for the subjet under consideration as is shown in Fig. 12. 

Here 

1. the primary quark emits the hard gluon with energy E in the direction 
of the angular interval 5, but not necessarily hitting the window, 

2. the emitted gluon produces the jet of partons with parton splitting 
angles larger than the window size, 

3. among those partons there exists such a parton (subjet) with energy k 
which hits the window, 

4. all decay products of that parton subjet cover exactly the bin 5. 

This picture dictates the rules of calculation of the q-th correlator of the 
whole jet. One should average the q-th correlator of the subjet F^'^^ over all 
possible ways of its production i.e. convolute it with the inclusive spectra of 
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such partons D'^^^ in the whole jet and with the probability of creation of the 
jet asKp. Analytically, it is represented by 

F^'HEoS) ~ f° 'f^mE/E,) f f D'«F<«)(M). (58) 

where Eq is the primary energy, E is the jet energy, and k is the energy of the 
subjet hitting the window. Since the unnormalized moments increase with 
energy while the parton spectrum decreases, the product D^^\k)F^^\k5) 
has a maximum at some energy fcmax) and the integral over momenta may 
be calculated by the steepest descent method. Leaving aside the details of 
calculations (see [25]), we describe the general structure of the correlator for 
the fixed couphng constant 70 = (Ga^/Tr)^/^ = const 

^ m{5)-'<'^'''{5Y''\ (59) 

where the three factors represent the phase space volume, the energy spec- 
trum factor and the g-th power of the average multiplicity. To get the nor- 
malized factorial moment Fq{5) one should divide (59) by the q-th. power of 
that part of the mean multiplicity of the whole jet which appears inside the 
window 5 i.e. by the share of the total average multiplicity corresponding to 
the phase space volume Ml: 

An(5) An{n). (60) 

If the analysis has been done in the d-dimensional space, the phase space 
volume is proportional to 

An ~ S^, (61) 

where S corresponds to the minimal linear size of the rf-dimensional win- 
dow. The last statement stems from the singular behaviour of parton prop- 
agators in quantum chromo dynamics (see [25]). That is why the factorial 
moments may be represented as products of the purely kinematical factor 
depending on the dimension of the analysed space and of the dynamical 
factor which is not related to the dimension and defined by the coupling 
constant 

~ (5)-'^(^-i)((5)(^'-')^°/«. (62) 
At small angular windows 6 the intermittency indices are given by 

0, = d(q - 1) - ^-^70- (63) 
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This formula is only valid for moderately small bins when the condition 
asln{A/6) < 1 is fulfilled. For extremely small windows, one should take 
into account that the QCD coupling constant is running. Then the constant 
7o should be replaced by the effective value (7) which depends logarithmi- 
cally on the width of the window S. As a result (see [25]), numerical values of 
the intcrmittency indices for very small bins become noticeably smaller than 
in the fixed coupling regime, especially for the low-rank moments. Moreover, 
the simple power-law behaviour is modified by the logarithmic correction 
terms and the intermittency indices depend on the value of the interval cho- 
sen. The resulting curve of In Fq{6) as a function of — In 5 has two branches 
and qualitatively reminds those shown in Fig. 1. The rather steep linear in- 
crease at the moderately small bins with the slope (63) is replaced at smaller 
windows by much slower quasi-linear increase. It is easy to calculate the 
location of the transition point to another regime and to show that at higher 
values of q it shifts to smaller bin sizes in accordance with trends in Fig. 1. 
We have described the results of the double logarithmic approximation of 
QCD. Higher-order terms have been treated in [25]. They do not spoil the 
general conclusions providing the corrections of the order of 10 per cents. 

It is interesting to note that higher dimension analysis just adds an in- 
teger number to the trivial part of the intermittency index and does not 
change its non-trivial "anomalous" dimension. Thus the increase of inter- 
mittency indices in higher dimension is trivial and has nothing to do with jet 
dynamics but is a consequence of phase space factors (mostly in the normal- 
ization denominator). The important difference from BE-effect is that the 
numerator of the moments provides the non-trivial "anomalous" part of the 
intermittency index which is absent in BE-treatment. It is common in all 
dimensions. 

The above results may be restated in terms of fractals. The power-hke 
behaviour of factorial moments points out to fractal properties of particles 

(partons) distributions in the phase space. According to the general theory of 
fractals (see [1] and references therein), the intermittency indices are related 
to fractal (Renyi) dimensions D{q) by the formula 

<P, = iq-l){d-Diq)), (64) 

wherefrom one gets (see (63)): 

D{q) = = 70 + -. (65) 
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The first term corresponds to monofractal behaviour and is due to the av- 
erage multiplicity increase. The second one provides multifractal properties 
and is related to the descent of the energy spectrum as discussed above. It 
is clearly seen that the fractality in quantum chromodynamics has a purely 
dynamical origin {D{q) ~ 70) related to the cascade nature of the process 
while the kincmatical factor in (64) has an integer dimension. The attempts 
to relate the fractal properties in the momentum space to the fractal struc- 
ture of coUiding objects in the ordinary space were tried also ([27]). The 
fractality in momentum space can be also formulated as the fractal nature 
of the subsequent available phase space at each branching of the gluon jet 



Coming back to our problem we would like to stress that the angular 
variable is the most convenient one to analyse correlations originated by the 
jet-like structure. At first sight the opening angle of the jet seems even 
a more convenient (and "natural") variable. However, it is more difficult 
to incorporate this angle into above theoretical study than just the usual 
polar angle 9. It is easy to show that the relative angle of two partons 6*12 
is connected with their polar emission angles ^1,^2, and with their relative 
distances in pseudorapidity 7^12 and azimuthal angle (fu by the formula 



at ^1 ^ 1, <^ 1- Therefore the analysis in the relative angle variable 
corresponds to two-dimensional analysis in r]i2 and Lpi2 when the size of the 
two-dimensional interval depends on polar angles of emitted partons and 
should be larger at small polar angles. It would produce higher average 
multiplicities in the denominator of factorial moments and suppress their 
values at small relative angles what is observed in experiment [14]. 

It exemplifies our statement about "natural" variables for each mechanism 
responsible for correlations in niultiparticle production. 

4 Discussion and conclusions. 

The very first intermittency studies were aimed at a scaling law in behaviour 
of factorial moments with hope to find out new collective effects in high 
energy interactions. Later it was recognized that it is just a part of the day- 
to-day work on correlation properties of multiparticle production that does 



([28, 29]). 




(66) 
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not diminish the importance of the above problem but helps also to disen- 
tangle contributions of different known mechanisms to particle correlations. 
According to our present-day theoretical prejudices we can name at least four 
of them. At the initial stage the quark-gluon jets appear. If described by 
the perturbative QCD they should give rise to the (quasi) intermittent be- 
haviour of factorial moments as functions of the (pscudo)rapidity bin size. At 
a simplified level of the hypothesis on local parton-hadron duality it should 
be valid for final hadrons as well. However, the transition from partons to 
hadrons could be not as simple as that, and it is sometimes considered as a 
phase transition imposing its own critical indices. In addition to it, the final 
stage interactions giving rise to known resonant states, surely, play impor- 
tant role. The final stage of creating hadrons asks for symmetry properties 
such as Bose-Einstein symmetrization to be respected. If one is interested in 
looking for new non-traditional sources of correlations (stochasticity, insta- 
bilities etc.) one should, first, to show that they contribute to correlations 
differently compared to considered "traditional" sources. 

As we mentioned above, each of them should be described in its own 
natural variable connected to its characteristic scale. We tried to show 
that resonances bring with them the mass scale of squared 4-momenta, BE- 
symmetrization is better revealed in 3-dimensional momentum analysis, while 
the jet-like structure asks for angular (or (pseudo)rapidity) variable with a 
scale determined by the corresponding "length" of the shower related to the 
(running) coupling constant. 

When analysed in "unnatural" variables, these mechanisms can produce 
the dependences which are not typical for them and mask (or imitate) some 
other effects due to the determinant of the transformation. It happens, 
for example, with BE-contribution when looked in Q^-variable. It becomes 
strongly peaked at small and imitates power-like law. One should not 
attempt to fit it by "traditional" Gaussian dependence when looked in that 
variable as well as one should not claim that it produces intermittency in a 
strict sense. 

Unfortunately, that example demonstrates also, that the contribution of 
some mechanism in its "unnatural" variable is not necessarily smeared out 
but, on the contrary, can produce rather steep dependence provoking mis- 
leading conclusions. 

The choice of (pseudo)rapidity as a "natural" variable in the original 
paper [2] was just related to traditions in theoretical approach and to exper- 
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imental facilities. It appeared to be natural for jets but not for resonances 
and B&correlations which seem to be smeared out in that variable. Sepa- 
rate study of the quantitative contributions of different mechanisms plotted 
as functions of the same variable asks for Monte-Carlo calculations. However, 
main qualitative ingredients are seen from above analytical approaches. 

In our opinion, the strong "kinematical" phase space dependence provided 
by the denominator of the normalized moments spoils the analysis introduc- 
ing strongly increasing (at small 5) factor depending on the dimension of 
the analysed bins. To unify the intermittency indices it looks appropriate 
to leave just (n) in the denominator of moments that cancels all kinematical 
factors at the expense of introducing energy dependence. 

To conclude, we have shown that different effects are better displayed if 
their "natural" variables are chosen. Some proposals are discussed above, but 
that study is just at the very initial stage and we call for further elaborated 
criteria beside those considered in the present paper. What concerns the title 
of the paper, we can say that Bose-Einstein correlations do contribute to 
intermittency in a wide sense. However, intermittency in its initial meaning 
exists even for unlike-charged particles and should be ascribed (probably) to 
jet-like parton cascades but not to Bose-Einstein correlations. 
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Figure captions 



Factorial moments of order q = 2,3,4 for the all charged 

sample, and the restriction to the positive-only and 

negative-only samples in ir'^p and K'^p interactions [9]. 

The rise of the factorial moments and density integrals 

with decreasing Sr) in pp collisions [10]. 

Log-log plot of the second factorial moment F2 in one 

dimension (a, b) and three dimensions (c, d) for unlike (a, c) 

and negative (6, d) charges. The full dots show the data, 

the open circles are Lund predictions [12]. 

Comparison of a), b) conventional factorial moments F2 

and F3 in three dimensions to c), d) F^^^ and 

F^^^ obtained from density integral method in GHP-topology. 

Sohd lines in a) and c) correspond to fits according to (25) [9] . 

The InF^^^^ as function of - ln((5VlGeV2). 

Note that the abscissa value 0.65 corresponds to the peak of 

p-meson and 1.77 is the value corresponding to the mass [9]. 

Log-log plot of a) F2{Q'^) for ( — ) and ( H — ) pairs and 

F3(g2) for ( ) and (ccc) triplets vs l/Q^ [15]. 

The rise of the density integrals with decreasing Q"^ in pp coUisions [10]. 
The normalized two-body correlation function for 
different charge combinations a) as a function of 5r], b) 
as a function of [10]. 

A comparison of pp collider data (UAl) with e+e" at the Z° 
pole (DELPHI) for the second order correlation function [17] . 

Log-log plot of K^iQ"^) for ( ) and (ccc) triplets a) from 

data and b) from the Lund Monte Carlo program including 
resonances but without BE correlations [15]. 
\n.Kq[Q'^) as a function of — InQ^ 

a) for like-charged and b) for unlike-charged particle combinations, 
compared to the expectations from FRITIOF model [18]. 
The subjet hitting the window 5 = 9 originates from 
a parton which appeared in the evolution of the primary 
gluon emitted by a parent quark. 
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